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Abstract 

The per-event yield of the highest transverse momentum charged particle and 
charged-particle jet, integrated above a given threshold starting at p^'^ = 0.8 
and 1 GeV, respectively, is studied in pp collisions at v^s = 8 TeV. The particles and 
the jets are measured in the pseudorapidity ranges \p\ < 2.4 and 1.9, respectively. 
The data are sensitive to the momentum scale at which parton densities saturate in 
the proton, to multiple partonic interactions, and other key aspects of the transition 
between the soft and hard QCD regimes in hadronic collisions. 
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1 Introduction 


The production of jets with large transverse momenta px ^ Aqcd ~ 0.2 GeV in high-energy 
proton-proton (pp) collisions originates from the scattering of partons, a process described by 
perturbative quantum chromodynamics (pQCD), through the convolution of the parton-parton 
cross section with the density of partons inside the protons. Jet production in pp collisions at 
the LHC, at transverse momenta pj > 20 GeV and in the pseudorapidity range | p | < 3, is well 
described by next-to-leading-order pQCD calculations [l]-EI- However, most of the final-state 
hadrons produced in pp collisions arise from the hadronisation of quarks and gluons scattered 
through "semi-hard" interactions with exchanged momenta of 0(1-3 GeV). At such low values 
of px, the theoretical partonic cross section, dtr /dpj oc iXg{pj)/pj, where is the strong cou¬ 
pling, becomes very large, and the integrated cross section (7(p™''^) = J^mm dpjda/dpj exceeds 

the total inelastic pp cross section, cpnei. At = 8 TeV, where cpnei ~ 70 mb |S, this occurs at 
p^'*^ values of 0(3GeV), much larger than the QCD scale, Aqcd/ at which the strong coupling 
diverges [Sjlll. 

Model calculations of hadronic collisions often regulate such an infrared divergence through 
an effective parameter connected to the confinement scale of hadrons fZl, such that the leading 
particle or leading jet production cross sections do not exceed the value of cpngi. Gontrary to the 
inclusive particle or jet production cross sections, the leading particle or leading jet production 
cross sections must indeed approach the total inelastic cross section because only one particle 
or one jet, the one with highest px in this case, is considered per event. In addition, at small px, 
the parton densities are probed in a region where parton recombination, i.e. saturation (see e.g. 
Ref. |8]), may occur. 

Reference l9l proposes that the jet cross section integrated over px > pl^ can be used as a 
probe of the transition from the perturbative (p!p™ » Aqcd) to the non-perturbative region 
(p^iri Aqcd)- According to Ref. El, this transition should also be visible for cross sections 
defined in restricted ranges of pseudorapidity. 

The results presented in this paper are based on measurements of single charged particles and 
jets reconstructed from charged particles alone. The advantage of jets is that they include more 
particles originating from the outgoing partons, while single charged hadrons carry only a 
fraction of the parent parton momentum. On the other hand, jets are sensitive to the underlying 
event (UE) activity, consisting of particles originating from multiple partonic interactions (MPI) 
and initial and final state radiation, while single leading tracks are not. The measurements 
based on leading particles and leading jets are therefore complementary. Throughout the text, 
the term "track-jets" refers to detector-level jets, reconstructed from charged-particle tracks 
observed in the detector, while "charged-particle jets" or just "jets" denote corrected, stable- 
particle level jets, consisting of stable charged particles from the final state. 

In this paper, the yields, r(p™™), for pp collisions with a leading charged particle or a leading 
jet are measured as a function of a minimum transverse momentum, 
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where Nevt is the number of selected events with a leading charged particle with pj > 0.4 GeV 
and \ri\ < 2.4, and N is the number of events with a leading charged particle or a leading jet 
with transverse momentum within |?/| < 2.4 or 1.9, respectively. 
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3 Experimental analysis 


2 Phenomenological models 

The measured distributions are compared to the predictions of different hadronic interaction 
models whose tunable parameters (mostly cormected to non-perturbative and semi-hard QCD 
phenomena) are obtained from comparisons to LHC data such as those on UE activity, inclusive 
multiparticle production and diffraction. 

The PYTHIA 6 I TOl and 8 ITTI event generators tame the low-px behaviour of the leading- 
order pQCD 2—>2 cross sections with a phenomenological factor HHHH [^':s(Pto + Pt)/^s(Pt)] 
[pf / (Pto + Pt)^]' where Pt,o is a (tunable) infrared regulator that runs with centre-of-mass en¬ 
ergy The tunes 4C fl^ , CUET llT^ . and MONASH IT^ are used, featuring different choices of 
the Pt,o cutoff, proton transverse profile, and/or parton distribution functions. 

The HERWIG++ liTSl Monte Carlo (MC) includes a hard (pQCD 2—>2 interactions) KT^ and a 
soft (non-perturbative) component for multiple interactions IITtI . The soft part is parametrised 
phenomenologically as dcr /dpj = Ae^^P^. The transition scale between the hard and the soft 
regions is set by the parameter px, 0 / obtained from fits to MPI and UE data, as well as to the 
effective cross section for double-parton scatterings. The parameters A and /3 are fixed by the 
requirements that the transverse momentum distribution be continuous at the matching scale 
px,0/ that the model reproduces the measured total cross section. HERWIG++ with tune 
CUETHSl [? ] is used to compared the measurement at reconstruction level. HERWIG++ is not 
used in the comparison of the final results, as it does not contain any diffractive component. 

The other two models, QGSJET-II [TSI and EPOS HU |20l, are based on Regge-Gribov effective 
field theory 11211 , which allows for a consistent treatment of soft and hard scattering processes in 
terms of the same degrees of freedom (reggeons and pomerons), based on unitarity cuts of the 
corresponding elastic scattering diagrams. Perturbative parton-parton processes are obtained 
via "cut (hard) pomeron" diagrams, and multi-scattering phenomena (saturation, MPI) are 
implemented through various procedures |[22ll . The two models differ in their approximations 
for the collision configurations (with exact energy sharing imposed in the case of EPOS), and 
the treatment of diffractive and perturbative contributions (the effective soft-hard transition 
occurs at px,o ~ 1.6 GeV for QGSJET-II and at px,o ~ 2GeV for EPOS). Pinally, in contrast to 
other MGs, EPOS includes also collective expansion effects in the final state that boost the final 
px distribution of the produced hadrons. It is worth to highlight that, for all MG models, the 
(centre-of-mass energy dependent) px,o cutoff plays a very similar role to the "saturation scale" 
(Qsat)/ which controls the onset of gluon fusion effects in the parton densities Ii23l . 


3 Experimental analysis 

The central feature of the GMS apparatus is a superconducting solenoid of 6 m internal di¬ 
ameter, providing a magnetic field of 3.8 T. A silicon pixel and strip tracker, a lead tungstate 
crystal electromagnetic calorimeter, and a brass and scintillator sampling hadron calorimeter 
are located within the volume of the solenoid. 

The inner silicon tracker measures charged-particle trajectories ("tracks" in the following) within 
the pseudorapidity range | p | < 2.5. It provides an impact parameter resolution of about 100 pm 
and a px resolution of about 0.7% for 1 GeV tracks at p = 0 112411 . A more detailed description of 
the GMS detector, together with definitions of the coordinate system and kinematic variables 
can be found in Ref. [25 |. 

The data analysed in this study were collected during a dedicated proton-proton run with an 
integrated luminosity of 45 pb^^ at a centre-of-mass energy of \/s = 8 TeV. This run has a low 
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instantaneous luminosity and a low probability (~2%) of multiple pp interactions occurring in 
the same bunch crossing (pileup). Pileup events are rejected by requiring exactly one vertex, 
following the method described in Ref. 12^ . 

Minimum bias events were selected online with the TOTEM T2 telescopes l27l that are placed 
symmetrically at about 14 m on both sides from the interaction point (IP). Single tracks are 
reconstructed in these telescopes with almost 100% efficiency for pj > 20 MeV/c, but because 
of multiple scattering and the effect of the magnetic field, tracks can be identified as coming 
from the IP with an efficiency that increases as a function of pj and is greater than 80% for 
Pj > 40MeV/c f28l . The minimum bias trigger, defined by the requirement of the presence of 
at least one track candidate in either of the T2 detectors [29l, has an efficiency close to 100% Il26l 
for events where a charged particle is produced within the T2 acceptance. According to the 
PYTHIA 8 and QGSJEtII- 04 |18l generators, about 91-96% of the total inelastic cross section at 
\/s = 8 TeV is seen by T2 JU, with the uncertainty coming mainly from low mass diffractive 
events. The present analysis follows the procedure described in Ref. |[26ll . where more details 
are given on the trigger, data selection, and correction procedures. 

Corrections for the contribution of background events triggered by T2 but without a charged 
primary particle in the T2 acceptance are estimated with simulated events from PYTHIA 8 
and EPOS. These models were found to enclose the measured pseudorapidity distributions 
of charged particles in the forward region |[26l . The average corrections for the two models 
vary from 4% and 1% at ~ 1 GeV to 7% and 5% at ~ 45 GeV, for the track and track- 
jet analysis, respectively. The deviation of PYTHIA 8 and EPOS from the average correction is 
taken as an estimate of the systematic uncertainty related to the T2 trigger efficiency; it is less 
than 0.7% for the leading track measurement and varies between 0.1 and 1.0% for the leading 
track-jet measurement Il26l . 

Events are selected offline by requiring the presence of a leading track in the region |p| < 2.4 
with pt > 0.4 GeV. These events are used to normalise the integrated distributions in both the 
leading track and the track-jet measurements. Track-jets are reconstructed offline from tracks 
with Pt > 0.1 GeV and |p| < 2.4, clustered by using the anti-kj algorithm 1301 - 13211 with a 
distance parameter of 0.5. The track-jet momentum is determined from the sum of all track 
momenta in the track-jet. The pseudorapidity restriction <1.9 assures that the track-jet is 
contained within the tracker acceptance. 

Detailed MG simulations of the GMS and T2 detectors are based on Geant4 1331 . Simulated 
events are processed and reconstructed in the same manner as collision data. Eor the correction 
of detector effects, as well as for comparison with models, both the PYTHIA 6 ITOl (version 
6.426) event generator with tune Z2* 1341 and the PYTHIA 8 (version 8.153) generator with tune 
4G are used. The final correction is obtained by averaging those from the two generators. 

The data are corrected to the stable-particle level, which is defined to include primary charged 
particles with lifetimes of cr > 1 cm, either directly produced in the pp collisions or from de¬ 
cays of particles with shorter lifetimes. According to this definition. Kg and A hadrons are con¬ 
sidered stable. Generated events are selected at the stable-particle level if at least one charged 
particle with pj > 40 MeV is present within the range 5.3 < \rj\ < 6.5, and at least one charged 
particle with pj > 0.4 GeV is found within |?/| < 2.4. In each event, the highest-px charged 
particle within |p| <2.4 and px > 0.8GeV is selected as the leading particle. Gharged particles 
are clustered into jets by using the anti-kx algorithm with a distance parameter of 0.5 with no 
restriction on px or p. The leading charged-particle jet is then defined as the charged-particle 
jet with the highest px above 1 GeV and < 1.9. 



4 


4 Results 


The average systematic uncertainty in the track reconstruction efficiency is taken to be 3.9% Il35l . 
Its effect is studied by randomly rejecting 3.9% of the tracks and then repeating the analysis. In 
the jet analysis, for tracks with low pj, the rejection probability is taken as 15% for px < 1 GeV. 
However, since the measurement is integrated over pj, it is nearly insensitive to even such large 
values of the rejection probability. The resulting uncertainty varies between 0.4% and 3.7% for 
the leading charged particle analysis and between 2% and 12% for the leading jet analysis. The 
larger uncertainties correspond to higher 


The Pj distribution of leading track-jets is unfolded to the stable-particle level by applying 
the iterative procedure Il36ll implemented in ROOUnfOLD Il37l in order to correct for the jet 
reconstruction efficiency and for migrations in jet pj. Thanks to the good pj resolution of the 
reconstructed tracks a simple correction for the track-finding efficiency is found to be sufficient 
for obtaining the pj distribution of leading charged particles. The PYTHIA 6 and pythia 8 MC 
models are used to generate the response matrices and efficiency corrections, and the average 
correction from the two generators is used to obtain the pj distributions at stable-particle level. 
The corrections vary between 5% and 10% at px ~ 1 GeV, to 10% and 40% at pj ~ 45 GeV, for 
the charged particle and the jet measurements, respectively. The deviation from the average 
is taken as an estimate of the systematic uncertainty related to the correction procedure. This 
uncertainty varies from 0.6 to 3% for the leading charged particle analysis, and from 2 to 10% 
for the leading jet analysis, depending on 

The systematic uncertainties are summarised in Table 

Table 1: The systematic uncertainties for the leading charged particle (0.8 < p^^^ < 50 GeV) 
and leading jet (1 < p!p™ < 50 GeV) measurements. 


Source 

Uncertainty (%) 

Leading charged particle 

Leading jet 

T2 trigger efficiency 

0.7 

O.l-l.O 

Tracking efficiency 

0.4-3.7 

2-12 

Gorrection procedure 

0.6-3.0 

2.0-10 

Total 

0.7-4.6 

2.5-16 


The per-event yields, defined in Eq. Q, are obtained experimentally as 


ripr) 


1 

Nevt 


E 



( 2 ) 


where Nevt is the number of events with a leading charged particle within \t]\ < 2.4 and with 
px > 0.4 GeV, Ap is the bin width, and AN is the number of events with a leading charged 
particle or leading jet in the bin. 


4 Results 

Figure shows the integrated distributions for the leading charged particle and leading jet 
events for p^^"^ >0.8 and 1 GeV, respectively. The distributions fall steeply at large transverse 
momenta, and by construction approach unity at small The turnover from a relatively flat 
to a steeply-falling distribution takes place between 1 and 10 GeV. However, the turnover point 
is different for the leading charged particles and the leading jet measurements. This reflects the 
fact that when particles are clustered into jets, more energy from additional particles is col¬ 
lected within the jet cone. In fact, when the jet cone size is reduced, the leading jet distribution 
approaches the leading charged particle distribution. 
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Figure 1: The integrated yield, of events with a leading charged particle within \ f]\ < 

2.4 (top) and with a leading jet within \ri\ < 1.9 (bottom), as a function of The data are 
compared to predictions from several PYTHIA 6 tunes (left) and various other event generators 
(right). The lower panels show the ratios of the MC and the data yields (MC/Data). The error 
bars indicate the statistical uncertainty and the red shaded area (only visible in the ratio plots) 
represents the systematic uncertainty. The predictions are scaled to the measured value of 
> 9.0 GeV) (top) and > 14.3 GeV) (bottom). The prediction from PYTHIA 6 with 

MPI off and no parton saturation is not shown in the MG/data ratio plot (left) because of the 
large disagreement with the data. 
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5 Summary 


For the comparison of the data to predictions of QCD MC generators the latter are rescaled 
to describe the high-pj‘*‘^ region. This rescaling is applied because the normalisation to the 
total visible cross section, which depends on the low-pj regularisation, affects the values of r 
also at where in fact theoretical predictions are more robust and agree better with 

the data. The exact choice of the normalisation point is arbitrary —> 9.0 GeV) for the 
leading charged particle, and > 14.3 GeV) for the leading jet—and the conclusions from 

this study are drawn from the shape of the distributions alone. The predictions at small 
thus give information on the modelling of the transition region from large to small 

In Fig. (left plots) the yields r(p^''^) as a function of p^^^ are compared to the predictions of 
the event generator pythia 6 with tunes Z2* and GUET, as well as with the default version 
of PYTHIA 6, both with and without MPI. Also shown is the impact of turning off the regu¬ 
larisation of the cross section, labeled "pythia 6 (default, MPI off, no sat)". At low p!p™, the 
distribution predicted by this latter model differs by more than one order of magnitude from 
predictions with the regularised cross section. 

In Fig. (right plots) the leading charged particle and leading jet data are compared with 
PYTHIA 8 with tunes 4G, GUET, and MONASH, HERWIG++ (version 2.7.0) with tune UE-EE-5G, 
EPOS (version 1.99) with LHG tune, and QGSJETII-04. 

The leading charged particle and leading jet cross sections are best described by EPOS, which 
deviates only by up to 10% from the data at very low p^''^ and reproduces the data well for 
p^™ > 4 GeV. The event generator HERWIG++ (UE-EE-5G tune) describes the leading jet cross 
sections fairly well, but does not reproduce the transition from large to small pj in the lead¬ 
ing charged particle cross section. The event generators PYTHIA 6 {Z2* and GUET tunes) and 
PYTHIA 8 (4G, GUET, and MONASH tunes) predict a somewhat different shape for the mea¬ 
sured distributions at small pj. 

The comparison of the MG predictions for MPI switched on and off indicates that the effect 
of MPI is small for leading charged particles, since the particle multiplicity plays only a minor 
role. However, when clustering particles into jets, the additional particles from MPI play a role, 
and a large difference is seen when such interactions are switched off in the simulation as in 
Fig. (bottom left); this brings PYTHIA 6 closer to the data at low p^*”. 

The predictions with MPI and saturation turned off (dashed curves in Fig.[^ left plots) exhibit 
a significant deviation from the data at small pj. 

In general, pythia and HERWIG++ describe the trend of the measured distributions but fail 
to reproduce the details in the 0(1-5 GeV) region, which calls for an improvement in their 
modelling of the transition from the non-perturbative to perturbative regime. 


5 Summary 

The integrated yields of events with a leading charged particle or a leading charged-particle 
jet with px above a given p^'*^ threshold, starting at p^™ = 0.8 and 1 GeV, respectively, have 
been measured in pp collisions at = 8 TeV in a data sample corresponding to an integrated 
luminosity of 45 pb^^. The particles and jets are measured in the pseudorapidity ranges |p | < 
2.4 and 1.9, respectively. 

The yields are found to be relatively flat in the p^^"^ region around 1 GeV—where the fixed- 
order perturbative parton-parton cross section diverges in the absence of any mechanism that 
saturates or unitarises the pQGD scattering—followed by a steep decrease for p^™ > 10 GeV. 
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The flattening behaviour observed at very low p™ is best described by EPOS, which deviates 
by at most 10% from the data. The comparison of the data with different phenomenological 
predictions of hadronic interaction models may help to improve the description of the transi¬ 
tion between the perturbative and non-perturbative QCD regimes, which is dominated by the 
effects of parton density saturation and multiple partonic interactions. 
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